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ABSTRACT 

We report Chandra ACIS-S3 x-ray imaging and spectroscopy of the R Aquarii bi- 
nary system that show a spatially resolved two-sided jet and an unresolved central 
source. This is the first published report of such an x-ray jet seen in an evolved 
stellar system comprised of -2-3 M© (also see Kellogg et al 2000). At E < 1 
keV, the x-ray jet extends both to the northeast and southwest relative to the cen- 
tral binary system. At 1 < E < 7. 1 keV, R Aqr is a point-like source centered on 
the star system. While both 3.5-cm radio continuum emission and x-ray emission 
appear coincident in projection and have maximum intensities at -7.5" NE of the 
central binary system, the next strongest x-ray component is located -30" SW of 
the central binary system and has no radio continuum counterpart. The x-ray jets 
are likely shock heated in the recent past, and are not in thermal equilibrium. The 
strongest SW x-ray jet component may have been shocked recently since there is 
no relic radio emission as expected from an older shock. At the position of the 
central binary, we detect x-ray emission below 1.6 keV consistent with blackbody 
emission at T - 2 x 10^ K. There is also a prominent 6.4 keV feature, a possible 
fluorescence or coUisionally excited Fe Ka line from an accretion disk or from 
the wind of the giant star. For this excitation to occur, there must be an unseen 
hard source of x-rays or particles in the immediate vicinity of the hot star. Such a 
source would be hidden from view by the surrounding edge-on accretion disk. 
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1. INTRODUCTION 



R Aqr is a symbiotic stellar system: a mass losing -1-2 M© Mira-like long period 
variable (LPV) with a 387 day period and a ~1 M© hot companion with many of 
its features explained by an accretion disk model. The companion is believed to 
give rise to the non-relativistic jet seen at UV, optical, and radio wavelengths. 
Since no UV continuum is observed, the accretion disk is believed to be edge-on 
(Kafatos, Michalitsianos, & HoUis 1986). The binary system orbit has been char- 
acterized as highly inclined to the line of sight (i ~ 70^) with large eccentricity (8 
~ 0.8), a semi-major axis of ~2.6xl0l4 cm, an orbital period of -44 yr, and at a 
distance of 200 pc (see HoUis, Pedelty, & Lyon 1997, and references therein). 
Hence, the northeast (NE)-southwest (SW) oriented jet is probably refueled epi- 
sodically with increased activity at periastron since the hot companion passes 
through the LPV's outer envelope. It is further posited that the binary system was 
at apastron circa 1996 (HoUis et al 1997a), so neither star should influence the 
other for many years near this epoch. 

From the early 1980s to the present epoch, the dominant NE radio jet has under- 
gone NE motion outward from the central HII region and counterclockwise rota- 
tion i.e. precession about the central HII region (HoUis, Pedelty, & Kafatos 
1997b). UV, optical, and radio observations of the NE jet are described as result- 
ing from shock excitation (e.g., HoUis et al. 1991; HoUis, Dorband, & Yusef- 
Zadeh 1992; HoUis et al 1997b) as the jet motion impacts the considerable cir- 
cumstellar material from LPV mass loss. For example, the thermal 
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bremsstrahlung radiation seen in the radio regime is suggested as a relic of a pass- 
ing shock front. Shock modeling of the R Aqr NE radio and UV jets suggests that 
temperatures should be on the order of that required for x-ray emission (HoUis et 
al 1997b). The R Aqr system was first reported as an unresolved Einstein x-ray 
source by Jura & Helfand (1984). It was later observed with EXOSAT (Viotti et 
al. 1987) and ROSAT (Hunsch et al 1998), so we planned a Chandra observation 
to determine the spatial and spectral structure and nature of the R Aqr x-ray emis- 
sion. 

2. OBSERVATIONS AND RESULTS 

2. 1 . Chandra Imagery and Spectra 

Chandra observed R Aqr for 24.5 ks on 2000 September 10 (JD 2451797.7 start 
time) with the Advanced CCD Imaging Spectrometer (ACIS)-S3 back-illuminated 
chip, which records photon energy for imaging and spectroscopy. The ACIS-S3 
pixel size is 0.492" (CXC 2000). The Chandra ACIS-S3 on-axis point-spread 
function in the 0.3 to 1.0 keV range has 50% and 80% encircled energy radii of 
-0.418" and -0.685", respectively (CXC 2000). A spatial image was constructed 
from the standard pipeline processed data. X-ray events were extracted from a 
0.38 arcmin^ region (5668 pixels) encompassing the jets and central source, 
yielding 1540 total counts in 22,717 s, or 0.068 counts s'l. 

Figure 1 is a contour image of smoothed Chandra data showing a resolved NE jet, 
an unresolved central source, and a resolved multi-component SW jet. The Figure 
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1 image was shifted in position by -0.4" in right ascension and -0.6" in declina- 
tion to align the peaks of the unresolved source in the x-ray and radio images (see 
§2.2 below). The Chandra absolute position uncertainty has a 1-a error circle ra- 
dius of 0.6" (CXC 2000) while the radio positions are accurate to 0.1", one tenth 
of the restoring beam size. 

Spectra of the NE jet, the central source, and the SW jet were extracted using 
CIAO Chandra analysis routines and XSPEC (Amaud 1996). Figure 2 shows the 
extraction regions, spectra and fits. In spite of small numbers of counts, the spec- 
tra required complex models for reasonable fits. A nonequilibrium coUisional ex- 
citation model was used for the jets (Borkowski et al 2001) with Anders & 
Grevesse (1989) (hereafter AG) initial abundances, with varying abundances for 
the SW jet (see Table 1). Data were binned into 15 counts ch"^ for the NE and 
SW jets and 114 eV channels for the central region. We fitted the interstellar col- 
umn density (the XSPEC WABS model). Results were consistent with radio ob- 
servations of 1.84 X 10^° cm"^ (Stark et al 1992). The NE jet shows two peaks, at 
-400 and 550 eV, and has few counts above 1 keV. The Rajmiond & Smith 
(1977) model gave a temperature in the 2 x 10^ K range with <1% confidence. It 
gave only a peak at the O VII line complex, -570 eV, increasing only to 6% con- 
fidence with variable abundances (the calculated N peak in that case was too 
broad to fit the data well, due to a complex of lines present in the model). There- 
fore, we used the NEI coUisional excitation XSPEC model (Borkowski et al 
2001). This gave an acceptable fit, including a second peak at -430 eV. The SW 
jet spectrum has only one peak at -550 eV. A Raymond-Smith fit gave a calcu- 
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lated peak at a higher energy than observed and only 7% confidence, not im- 
proved by varying the abundances from the AG values. The XSPEC mekal mod- 
els gave similar results. A good fit was obtained using the NEI model only with 
relative O abundances >20. At relative O abundances < 20, the calculated O peak 
was at a higher energy than the data peak. More complex coUisional excitation 
models having more free parameters did not improve the fit. 

In Figure 2, we also show a fit to a blackbody for the central source at low ener- 
gies, and nine gaussians at expected K line energies of Fe and the lighter even Z 
elements, down to 1.6 keV, where the blackbody model dominates. We quote the 
blackbody fit parameters in Table 1, although we could not reject equilibrium op- 
tically thin hot gas models at the 50% confidence level. The hard central source 
shows a feature containing 16 counts with a flux of 4.9 ± 1.7 x lO'^^ ergs cm'^ s' 
1 at 6.41 ± 0.04 keV. This could be Fe I (i.e., Ka), or up to Fe IV, based on the 
calculations tabulated in Kaastra and Mewe (1993), or up to Fe X based on the 
measurements of Decaux et al (1997). At very low statistical significance, there 
are also suggestions of emission from Fe KP at 7.1 keV, 3 counts, and emission 
from Ka lines of Ca, V, and Cr: 7 counts at 3.7, 6 counts at 5.0, and 4 counts at 
5.4 keV, respectively. The expected number of background counts in each chan- 
nel is a few x 10" ; there are essentially no background counts in this energy 
range. This was checked by displacing the central source region about 1"; it con- 
tained no counts. Our data are very likely real x-ray emission events from the cen- 
tral source. 
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We observed a flux of 1.0 ± 0.2 x 10'^^ erg cm"^ s'^ at E < 2 keV from R Aqr. 
Hunsch et al (1998) observed 5x greater flux in that energy range during the 
ROSAT all sky survey, epoch 1990-1991. We also observed an additional 1.0 ± 
0.25 X 10"^^ erg cm"^ s'^ at E > 2 keV, not visible to ROSAT. Viotti et al (1987) in- 

11 2 1 

f erred ~2 x 10" erg cm" s" by estimating the spectrum of the source. They also 
interpreted the Einstein report as an upper limit that corresponds to ~ 2.5 x 10" 
erg cm"^ s"\ We conclude that R Aqr is several times fainter in our data than it 
was in the 1980s and early 1990s. 

2.2. VLA Imagery 

R Aqr was observed on 1999 November 1, 4, and 7 with the NRAO^ Very Large 
Array in the B configuration. Two sets of receivers were combined to observe a 
100 MHz band centered at a frequency of 8460. 1 MHz (-3.5 cm wavelength). A 
total of ~2 million visibilities were accumulated in ~8 hours of on-source integra- 
tion. AIPS software was used for standard calibration and imaging. The resulting 
uniformly weighted image is thermal noise limited (la noise level of 8 |lJy) and 
has a restoring beam of 1.00" by 0.69" with a major-axis position angle of -3.9°. 
The NE radio jet is shown as contours in Figure 3 overlaid on a color image of the 
NE x-ray jet, and both jets are seen in projection to be spatially coincident. 



^ The National Radio Astronomy Observatory is operated by Associated Universities, Inc., under 
cooperative agreement with the National Science Foundation. 



10 



3. DISCUSSION 

The R Aqr symbiotic star system is strikingly different at E < 1 keV than from 1-7 
keV. At the lower energy, the jets dominate. Above 1 keV, R Aqr looks point-like 
at the Chandra resolution, centered on the star system. The x-ray jets are another 
manifestation of directed flow seen in the radio, optical, and UV. This is solid 
evidence for cylindrical accretion geometry, i.e. a disk. HoUis and Koupelis 
(2000) developed a parcel ejection model for the R Aqr jet emitted along the rota- 
tion axis of the massive rotating star-accretion disk system that anchors a mag- 
netic field. A natural extension of the model suggests that the x-ray emission 
from the jet could be the result of shock emission as the jet parcel plows into the 
ambient circumstellar medium. The preservation of a well-defined direction for 
the jets, with evidence for precession, argues against simple spherical accretion. 
In Figure 1, the strongest parcel in the SW jet seems to have contour compression 
on its SW edge as would result from SW directed motion. The NE jet component 
has contour compression on its eastern edge that may be due to counter clockwise 
motion about the central source reported previously for the NE radio jet (HoUis et 
al 1997b). In Figure 3, the NE x-ray jet as seen in projection lies along the same 
trajectory as the NE radio jet, supporting a counter clockwise rotation scenario for 
the NE x-ray jet. The 430 eV peak in the NE jet and the need for the NEI model to 
fit the spectrum suggests that this feature was formed recently, with a product of 
age and ambient density of 4.3 x 10 s cm" . If n ~1 cm" as for the typical inter- 
stellar medium, the age is -15 years. This is comparable to the observed timescale 
for changes in the radio emission, and in x-ray, as inferred from comparing our 
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data with the earlier x-ray satellite observations. For greater density, as might be 
expected near the R Aqr Mira, the age may be less. The SW x-ray jet is -30" from 
the central source compared with -7.5" for the NE jet, suggesting that it was 
ejected earlier. It has a different spectrum, with no N peak. We could only fit a 
model with a very high O abundance. Either there is a velocity shift in the O peak, 
or there is a different mix of the O VII lines at 0.569, 0.574 and 0.666 keV than 
our models can reproduce. The derived value of the ionization timescale is about 
half that of the NE jet, so given the expected variation in ISM density, the time 
since heating may be about the same. However, if the SW jet was shocked further 
in the past, we should expect to see relic radio emission as in the NE jet that is not 
detected. The shock temperatures from the two NEI fits to the jets seem disparate, 
but only differ by -1.5a, so the difference is not very significant. The central x- 
ray source has a counterpart in 3.5-cm radio continuum (c.f. Figure 3). The radio 
emission is a thermal HIT region at -10^ K while the central x-ray source low en- 
ergy spectrum in Figure 2 fits indicates a temperature of order 10^ K. Perhaps the 
x-ray emitting volume is interior to the Hll region, associated with a high energy 
density region at the hot star/accretion disk. The neutral or lightly ionized Fe at 
6.4 keV, without strong associated continuum is evidence for relatively cooler 
material excited by an unseen source of hard photons or charged particles. Such 
emission could be reflecting from the surface of the edge-on accretion disk (e.g., 
Frank, King, & Raine 1992, p. 203), or directly exciting the relatively cool mate- 
rial at the boundary of a surrounding Stromgren surface. 



IV. SUMMARY 
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We have imaged a two-sided x-ray jet and detected an unresolved central x-ray 
source in the R Aqr binary system with Chandra ACIS-S3. Non-equilibrium col- 
lisional ionization models were applied to the jets. Both jets have ionization time- 
scales of -10*^ s cm"^, corresponding to ages of order 10 y or less. A 6.4 keV fea- 
ture in the central source must be a de-excitation line of relatively cool gas, evi- 
dence of an unseen hard source near the hot star, hidden from view owing to the 
surrounding edge-on accretion disk. Comparison of the NE x-ray jet and central 
x-ray source with 3.5-cm radio continuum VLA imagery at the same resolution 
suggests that the x-ray jet may be a high temperature, low density region associ- 
ated with cooler post-shock UV and radio emitting regions. After more than one 
hundred years of intense study, R Aqr is still revealing new and fascinating fea- 
tures. 

We thank Drs. Eric Feigelson , Steve Drake, John Raymond and Brad Wargelin 
for useful discussions, and the referee for a thorough reading of the manuscript 
and suggestions for improvement. We acknowledge support from NASA grant 
GO0-1062A. E.K. received support from NASA contract NAS8-39073. 
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Table 1. Spectral fit parameters^ 



Jets: Non-equilibrium ionization model with interstellar absorption 



Quantity 


NE Jet 


SW Jet 


Counts m spectrum, E > 0.3 keV 


396 


196 


C abundance (Fraction of AG) 


i.U (iixea) 


< 0.5 


N abundance (Fraction of AG) 


1.0 (fixed) 


0. (fixed) 


O abundance (Fraction of AG) 


1.0 (fixed) 


2 (fixed) 


nH(10^"cm"^^ 


< 5.4 ± 1.5 


< 10 


kT (keV) 


1.66 ±0.99 


0.20 ± 0.06 


T(K) 


19.3 ± 11.5 


2.3 ±0.7 


X, ionization timescale (s cm' ) 


4.3 ±0.3 X 10*^ 


1.9 ±0.5 X 10^ 


Normalization (10"^ ph cm'^s'' keV') 


4.0 ± 1.0 X 10'^ 


1.5(0.8-7.1) X 10"^ 


Observed flux (10"'"* ergs cm"^ s"^) 


7.1 ± 1.8 


2.6(1.4-12.3) 


Source flux (nH= 0, 10"^"^ ergs cm"^ s"^) 


18.9 + 4.7 


3.1(1.7-14.7) 


Reduced %^ 


1.04 


1.19 


Confidence (%f 


41 


31 



Center soft source: Blackbody with interstellar absorption 



Quantity 


Center Source 


Counts in spectrum, E > 0.3 keV 


108 


nH(10^"cm-^) 


<6.2 


kT (keV) 


0.18 ±0.02 


T (lO^K) 


2.09 ± 0.23 


Normalization (10"^ ph cm"'^ s"^ keV"^^ 


0.009 ± 0.005 


Observed Flux (10'^^ ergs cm"^ s"\ 0.3-2.0 keV) 


0.68 ±0.38 


Source Flux, nH= (10"'^ ergs cm"^ s\ 0.3-2.0 keV) 


0.68 ±0.38 


Reduced %^ 


0.55 


Confidence (%) 


84 



^ AH errors are 1 a. 

^ XSPEC null hypothesis probabiUty. 
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FIGURE CAPTIONS 

Fig. 1. Contour map of the R Aqr 0.1-10 keV x-ray image showing the NE jet, 
central source, and SW jet. The original ACIS image has been convolved with a 
1" FWHM circular gaussian to facilitate smooth contours. The peak photon count 
in this image is 194.4 counts and the contour levels are 1.4, 2.8, 4.2, 5.6, 8.4, 11.2, 
16.8, 22.4, 33.6, 44.8, 89.6, and 179.2 counts. At E>1 keV, the jets are not visible; 
the dominant emission comes from a point-like source at the star system. 

Fig. 2. R Aqr x-ray spectrum extraction regions and spectral fits. The upper panel 
shows the raw image with source and background spectral extraction regions; 
NE2a and Bla for the NE jet, cntr and B2a for the central source, and SWa and 
B2a for the SW jet. The plots show the spectra of the NE jet, central source, and 
SW jet with best fit models. The NE jet data, binned to 15 counts per channel, 
has peaks from expected O VII emission at 0.574 keV and N VI at 0.43 keV. The 
histogram fit is a non-equilibrium ionization model with interstellar absorption. 
The central source data, binned in 1 14 eV width channels, is modeled by a black- 
body, a Fe Ka line at -6.4 keV, Kp at -7.0 keV and other discrete but weak lines 
at K energies of Cr, V, Ca, Ar, and S. The Fe line suggests a hidden hard source 
associated with an accretion disk. The SW x-ray jet data, binned to 15 counts per 
channel, has the expected O VII peak at 0.574 keV. The fit is to a non-equilibrium 
ionization model with high O abundance, (e.g., see Table 4 of HoUis et al 1991). 
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Fig. 3 - Color rendering of the Chandra ACIS image of the R Aqr NE jet and cen- 
tral HII region overlaid with a contour plot of the 3. 5 -cm radio continuum emis- 
sion, observed ten months earlier. The coordinate system is defined relative to the 
peaks of the x-ray and radio continuum emission (assumed coincident). The x-ray 
image has been convolved with a 1.00" x 0.69" elliptical gaussian that is the same 
size as the restoring beam of the radio image. The peak flux is 10.94 mJy beam'l 
and the contour levels are 0.02, 0.04, 0.06, 0.08, 0.12, 0.16, 0.24, 0.32, 0.48, 0.64, 
1.28, 2.56, 5.12, 10.24 mJy beam"!. 
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